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We demonstrate the implementation of a quantum al- 
gorithm for estimating the number of matching items in a 
search operation using a two qubit nuclear magnètic reso- 
nance (NMR) quantum computer. 



Quantum computers |ï],|| offer the tantalising prospect 
of solving computational problems which are intractable 
for classical computers. A variety of algorithms have 
been developed, most notably Shor's algorithm for fac- 
torising composite numbers in polynomial time |ï|,Q , and 
Grover's quantum search algorithm (^||. Until recently 
these algorithms were only of theoretical interest, as it 
proved extremely difhcult to build a quantum computer. 
In the last few years, however, there has been substantial 
progress |7]-|í| in the construction of small quantum com- 
puters based on nuclear magnètic resonance (NMR) stud- 
ics [|Ï0| of the nuclei of small molècules in solution. NMR 
quantum computers have been used to implement a va- 
riety of simple quantum algorithms, including Deutsch's 
algorithm |ïï| , [ï^| and Grover's algorithm fÏ3| - |ï^|. 

NMR quantum computers differ from other implemen- 
tations in one important way: there is not one single 
quantum computer, but rather a statistical ensemble of 
them. For this reason NMR quantum computers should 
be dcscribcd using dcnsity matriccs rather than the more 
usual ket notation. In some cases this ensemble nature 
is irrelevant: it is possible to prepare the system with 
an initial density matrix indistinguishable from that of 
a pure eigenstate (a pseudo-pure state), and as long as 
the result is another pseudo-pure state the behaviour of 
an ensemble quantum computer is identical to that of a 
conventional quantum computer. Some algorithms, how- 
ever, produce a superposition of states (relative to the 
natural NMR computational basis) as their final result, 
and in such cases the behaviour of an ensemble quantum 
computer will be quite different. 

An important example is Grover's algorithm whcn 
there is more than one matching item to be found | {Ï6| . 
Suppose a search is made over TV items amongst which 
there are k matching items. After 0(y/N/k) evaluations 
of Grover's search function the quantum search algorithm 
will produce an equally weighted superposition of the k 



matching items. With a conventional quantum computer 
this state allows any one of the k matching items to be 
determined at random, as a measurement will result in 
one of the states contributing to the superposition. With 
an ensemble quantum computer, however, different mem- 
bers of the ensemble will result in different states, and the 
final observed signal will be an average over the k match- 
ing vàlues. In general it will be difficult or impossible to 
deduce anything about individual matching items from 
this ensemble average, and so NMR quantum computers 
will not be capable of carrying out conventional Grover 
searches when more than one item matches the search 
criteria. 

An alternative approach to searching is to count the 
number of matching items found in some desired portion 
of the search space. Clearly a bisection search will then 
permit the first matching item, for example, to be lo- 
cated in approximately log 2 (iV) attempts. This is only a 
sensible strategy if some efficient algorithm for counting 
matches can be found. Fortunately this can be achieved 
by a simple modification of Grover's quantum search, ap- 
proximate quantum counting [p^|-p^[ . 

Suppose we have a function f(x) which maps n-bit bi- 
nary strings to a single output bit, so that f(x) = or 
1. In general there will be TV — 2™ possible input vàl- 
ues, with k vàlues for which f(x) — 1. Grover's quantum 
search j||||ï(J allows one of these k items to be found, 
while quantum counting 16 18) allows the value of k 
to be estimated. The counting algorithm can be con- 
sidered as a method for estimating an eigenvalue of the 
Grover iterate G — HUqH~ 1 Uj, which forms the basis of 
the searching algorithm (the operator H corresponds to 
the n-bit Hadamard transform, Uq maps |0) to — 10) and 
leaves the remaining basis states alone, and Uj maps \x) 

to (-l)^) +1 |x)). 

Starting from the state 1 000 ... 0) (000 . . . 0| apply the 
Hadamard operator H to obtain an equally weighted su- 
perposition of all basis states. For < k < N we write 



ff|Q0Q...0) = (|* 
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where |^+) and l^-) are eigenvectors of G with eigenval- 
ues e and cos(0fc) = 1 — 2k/N. For the two extreme 
cases, k = and k — N, H |000 ... 0) is itself an eigenvec- 



tor, and we can write yS 



+/ 



i2]QOQ...O), with 



eigenvalues given by the formulae above. 
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Eigenvalue estimation is most easily described by con- 
sidering a register which begins the calculation in an 
eigenvector of G, say \^+). An additional control qubit is 
needed which begins in the state (|0) + |1))/V^; this may 
be obtained from |0) by a Hadamard transform. The op- 
erator G is then applicd to the target register when the 
control bit is in state |1), that is, a controlled-G. The 
controlled-G produces the result 



v/2 



(|0)+e i *»|l))|*. 



or after r repetitions of the controlled-G 
Í=(|0)+e^|l»|* + >. 
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Applying a second Hadamard transform to the control 
qubit gives 
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(4) 



tracing out the target register and expanding the expo- 
nential terms gives for the final sate of the control qubit 



1 /l + cos(r0 fc ) isui{r(t> k ) 



2\ -isin(r<j) k ) 1 - cos (r </> fc ) 



(5) 



The same result is obtained if we replace 1*+) with \^-), 
except that the two off diagonal elements are negated. 
Thus the same diagonal elements are also obtained from 
any superposition or statistical mixture of the two, such 
as ií| 000... 0) (equation ([]). 

A variety of different ensemble measurements can be 
performed to characterise the final state of the control 
qubit, but the simplest approach is to measure the ex- 
pectation value of a z . This corresponds to determining 
the population difference between the |0)(0| and |1)(1| 
states, and is proportional to cos(r<pk)- Note that in this 
case ensemble quantum computers have an an advantage: 
with a single quantum computer it would be necessary 
to repeat the calculation several times in order to obtain 
a statistical estimate of cos(r</> fc ). 

<j>k can be estimated by varying r (the number of repe- 
titions of the controlled-G). Estimating <f>k with sufficicnt 
accuracy to determine k requires roughly yJk(N — k) ap- 
plications of G (ïl| , while a classical algorithm would re- 
quire N evaluations of /. It is also possible to estimate k 
to some desired accuracy: to obtain an estimate k with 
accuracy e, that is 



\k - k\ < ek 



(0) 



requires on the order of (l/e)y/N/(k + 1) applications of 
G JÏ7| , fL8|]2"o"| , while a classical algorithm requires about 
(l/e*)N/(k + 1) evaluations of /. 

A quantum circuit for implementing this algorithm on 
a two qubit NMR quantum computer is shown in fig- 
ure ^. This differs from the conventional circuit in two 



ways. Firstly pairs of Hadamard gates are replaced by 
an NMR pseudo-Hadamard gate (a 90° rotation) and its 
inverse |ïïj . Secondly the controlled-Hadamard gates in- 
side the controlled-G propagator have been replaced by 
uncontrolled gates; this is permitted as the intervening 
Uq gate has no effect when the control spin is in state 
|0). This circuit can be used to count the number of 
solutions to f(x) = 1 over a one bit search space, but 
similar circuits exist for larger search spaces. 
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FIG. 1. A quantum circuit for implementing quantum 
counting on a two qubit NMR quantum computer; the central 
sequence of gates, surrounded by brackets, is applied r times. 
The upper line corresponds to the control bit, while the lower 
line corresponds to the target bit. A similar circuit can be 
constructed for a larger search space by replacing the target 
bit by a register and replacing gates applied to the target 
by multi-bit versions. Gates marked h implement the NMR 
pseudo-Hadamard operation, while those marked h' 1 imple- 
ment the inverse operation. Controlled gates are marked by 
a circle and a vertical "control line" . 

This algorithm was implemented using our two-qubit 
NMR quantum computer O], which uses two l· l·l nuclei 
in a solution of the small molecule cytosine in D2O. All 
NMR experiments were carried out on a home-built spec- 
trometer at the Oxford Centre for Molecular Sciences, 
with a l· l·l operating frequeney of 500 MHz. The two spin- 
states of the l· l·l nuclei act as qubits, and it is necessary 
to address each spin individually. Previous experiments 
on this system |ÏÏ],[Ï4| have used soft pulses to achieve 
selective excitation, and errors in these pulses have re- 
sulted in significant distortions in observed spectra. For 
these experiments a different approach was adopted, us- 
ing non-selective hard pulses whenever possible. 

The l· l·l transmitter frequeney was set in the centre of 
the spectrum, so that the two spins have angular freqüèn- 
cies in the rotating frame of ±w/2. The Hamiltonian can 
then be written in product operator notation |23| as 



H = —I z — —S z + ttJis 2I z S z 



(7) 



where Jis is the spin-spin coupling constant, and weak 
coupling has been assumed (i.e., lo 3> Jis)- Using a com- 
bination of non-selective pulses and carefully chosen peri- 
ods of free evolution under Tí it is possible to implement 
many of the necessary gates without the use of selective 
pulses. For example the controlled-C/-^- gate, which im- 
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plements the function when /(0) = and /(l) = 1, can 
be constructed using the pulse sequence shown in figure]^. 



180 180 

S + £270 + 8 X 5 



FIG. 2. A pulse sequence implementing a controlled-C/^- 
gate using only hard pulses and periods of free precession. 
Pulse rotation angles (in degrees) are marked above each 
pulse, while pulse phases are marked within a pulse. Other 
periods correspond to free precession under the Hamiltonian 
Tí for the time indicated. These times are chosen such that 
45 + £270 = 1/(2 J/s) and e 2 70 = Stt/lü. 

Some gates, however, cannot be implemented without 
using selective pulses; for example the pseudo-Hadamard 
gates within the controlled-G should only be applied to 
the target spin. Fortunately it is possible to create se- 
lective pulses using only hard pulses and delays, and 
this process is particularly simple when only two spins 
are involved. For short periods of evolution under Tí 
the small spin-spin coupling term can be neglected, and 
Tí w (u>/2)(I z — S z ). Thus after a time £45 = ir/(2uS) the 
two spins will have undergone rotations of ±45° about 
their respective z-axes. This ±z-rotation can be con- 
verted to a ±j/-rotation by sandwiching the r period be- 
tween 90° and Q0°_ x pulses (a variant of the more tra- 
ditional composite z-pulse pi). Combining this with a 
45° pulse along the y-axis gives an overall 90° rotation 
for the first spin (/), but no nett rotation for the second 
spin (S), as shown in figure]^. 
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FIG. 3. The sequence of hard pulses and delays shown in 
(a) is equivalent to the single selective pulse (b); other se- 
lective pulses can be implemented in a similar fashion. Note 
that the small spin-spin coupling may be neglected during the 
short period £45. 

With minor variations this approach can be used to 

generate selective pulses along any desired axis, and 

which excite eithcr I or S as desired. These selective 

pulses can then be used to implement the remaining 

gates: for example a controlled-C/j-^- can be implemented 

using the circuit for controlled-[/-^- with a selective 180° 

... J 01 
pulse applied immediately before and after the other 

pulses. 

The circuit shown in figure ^ encodes the result of the 
calculation in the state of the control qubit. This state 



could be characterised in a variety of ways, of which 
the simplest is to measure the expectation value of a z 
for the spin. This cannot be achieved directly, as z- 
magnetisation is not a direct NMR observable, but an 
equivalent measurcmcnt can be easily made by exciting 
the spin with a 90° pulse and then observing the result - 
ing NMR spectrum. After appropriate phase correction 
the integrated intensity of the corresponding signal gives 
the desired result. The phase correction step requires a 
reference spectrum , but this is easily obtained by 

acquiring a spectrum with r = 0. 

Immediately prior to acquisition a short magnètic field 
gradient pulse was applied to destroy the homogeneity 
of the main field. This has the effect of dephasing (and 
thus rendering undetectable) all off-diagonal terms in the 
final density matrix Q] , with the exception of those cor- 
responding to zero quantum coherence fï(Í|j . The zero 
quantum terms can also be removed using the fact that 
they evolve at freqüències of ±lu under the Hamiltonian 
Tí. This zero quantum filter is easily combined with a 
Standard four-step CYCLOPS phase-cycle |üj, to reduce 
instrumental imperfections. 




FIG. 4. Experimental results from our NMR quantum com- 
puter for each of the four possible functions, /. The observed 
signal intensity is plotted as a function of r, the number of 
times the controlled-G operator is applied, and all intensities 
are normalised relative to the case of r — 0. The sòlid lines 
are exponentially damped cosinusoids with the theoretically 
predicted freqüències, and are plotted merely to guide the eye. 

The results of our NMR experiments are shown in fig- 
ure f|. Measurements were made for each of the four 
possible functions: / 00 (k = 0), / i (k = 1), fio (k = 1), 
and fií (k = 2). In each case the predicted signal is a 
cosinusoidal modulation of the signal intensity as a func- 
tion of r, the number of times the controlled-G is applied, 
where the frequency of the modulation, <pk, depends on 
k/N. For the two-qubit case, where N = 2, the behaviour 
is particularly simple, with modulation freqüències of 
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(k = 0), n/2 (k — 1), and ir (k = 2). In this case it is 
possible to determine k using just one experiment, with 
r = 1, but spectra were also acquired with larger vàlues 
of r, both to demonstrate the principle involved and to 
explore the build up of errors in the calculation. 

The experimental results do indeed show a cosinusoidal 
modulation as expected, but they deviate from the sim- 
ple predictions above in a number of ways. Firstly all 
the signals show a clear decay in signal intensity as r is 
increased, and this decay is most ràpid for /oi and fio 
(where fc = 1), and least ràpid for /n (where k = 2). 
The simplest explanation for this observation is decoher- 
ence: for large vàlues of r the total length of the pulse 
sequence is comparable to the spin-spin relaxation time, 
T2. Another likely cause is imperfections in the pulses 
applicd, in particular those arising from variations in the 
strength of the resonant RF field across the sample (B\ 
inhomogeneity) . Both effects are expected to be most 
severe when k — 1 , as these cases have complex Uj gates 
which take a long time to implement, and least severe 
when k — 2, in which case Uj is just the identity opera- 
tion. 

In addition to the main exponential decay other de- 
viations from the simple behaviour predicted by theory 
can be seen. These effects are clearest for /oo (k = 0), 
where alternate signal intensities are seen to lie alter- 
nately above and below the main curve. Such effects 
could in principle arise from many different causes, but 
numerical simulations indicate that the major cause is 
off-resonance effects. These occur because the applied 
RF field is not perfectly resonant with the NMR tran- 
sitions, but instead is applied a small distance (±o;/2) 
away. Thus the effect of the field (in the rotating frame) 
is not simply to cause a rotation around itself, but rather 
to cause a rotation around a tilted axis |u|. We are 
currently seeking ways to reduce the size of such effects. 

Despite these small errors the results are remarkably 
good, especially for the case of f\\. In this case the ex- 
periments have been repeated with much larger vàlues 
of r, and the cosinusoidal variation remains clearly visi- 
ble after 60 or more iterations (data not shown). Thus 
our NMR quantum computer is capable of demonstrating 
quantum algorithms involving several hundred quantum 
gates. 
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